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ABSTRACT 

Recent measurements of type la supernovae as well as other concordant ob- 
servations suggest that the expansion of our universe is accelerating. A dark 
energy component has usually been invoked as the most feasible mechanism for 
the acceleration. However, the effects arising from possible extra dimensions can 
mimic well the role of a dark energy through a modified Friedmann equation. In 
this work, we investigate some observational constraints on a scenario in which 
this modification is given by = ^y^(p + Cp""). We mainly focus our attention 
on the constraints from recent measurements of the dimensionless coordinate dis- 
tances to type la supernovae and Fanaroff- Riley type lib radio galaxies compiled 
by Daly and Djorgovski (2003) and the X-ray gas mass fractions in clusters of 
galaxies published by Allen et al. (2002,2003). We obtain the confidence region 
on the power index n of the modificative term and the density parameter Qrn 
of the universe from a combined analysis of these databases. It is found that 
n = 0.06lo;^g and fim = 0.301q;o25 the 95.4% confidence level, which is consis- 
tent within the errors with the standard ACDM model. These parameter ranges 
give a universe whose expansion swithes from deceleration to acceleration at a 
redshift between 0.52 to 0.73. 
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1. 



Introduction 



The Hubble expansion, the Cosmic Microwave Background Radiation(CMBR), the pri- 
mordial Big Bang Nucleosynthesis and the structure formation are the four pillars of the 
standard Big Bang cosmology. Recent years, it seems that all these cornerstoness combined 
to point out that the expansion of the universe is speeding up rather than slowing down (for 
a recent review see Peebles and Ratra 2003). The main evidence comes from the recent well 
known distance measurements of some distant Type la supernovae (Pcrlmutter et al. 1998, 
1999; Riess et al. 1998, 2001). Possible explanations for such an acceleration include: a cos- 
mological constant A (Weinberg 1989; Carroll et al. 1992; Krauss and Turner 1995; Ostriker 
and Steinhardt 1995; Chiba and Yoshii 1999), a decaying vacuum energy density or a time 
varying A-term (Ozer and Taha 1987; Vishwakarma 2001; Alcaniz and Maia 2003; Jain, Dev 
and Alcaniz 2003), an evolving scalar field (referred to by some as quintessence: Ratra and 
Peebles 1988; Caldwell et al. 1998; Wang and Lovelace 2001; Li, Hao and Liu 2002; Weller 
and Albrech 2002; Li et al. 2002a,b; Chen and Ratra 2003; Mukherjee et al. 2003), the 
phantom energy, in which the sum of the pressure and energy density is negative (Caldwell 
2002; Hao and Li 2003a,b; Dabrowski et al. 2003), the so-called "X-matter" (Turner and 
White 1997; Zhu 1998,2000; Waga and Miceh 1999; Podariu and Ratra 2001; Zhu, Fujimoto 
and Tatsumi 2001; Sereno 2002; Alcaniz, Lima and Cunha 2003; Lima, Cunha and Alcaniz 
2003), and the Chaplygin gas (Kamenshchik et al. 2001; Bento et al. 2002; Alam et al. 
2003; Alcaniz, Jain and Dev 2003; Dev, Alcaniz and Jain 2003; Silva and Bertolami 2003; 
Makler et al. 2003). All the above mechanisms for accelerating are obtained by introducing 
a new hypothetical energy component with negative pressure - the dark energy. 

On the hand, many models have appeared that make use of the very ideas of branes 
and extra dimensions to obtain an accelerating universe (Randall and Sundrum 1999a,b; 
Deffayet et al. 2002; Avelino and Martins 2002; Alcaniz, Jain and Dev 2002; Jain, Dev and 
Alcaniz 2002). The basic idea behind these braneworld cosmologies is that our observable 
universe might be a surface or a brane embedded in a higher dimensional bulk spacctime in 
which gravity could spread (Randall 2002). The bulk gravity see its own curvature term on 
the branc which accelerate the universe without dark energy. Here we are concerned with 
the cosmological model from the modified Fricdmann equation as follows 
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Preese and Lewis (2002) showed the above term proportional to (dubbed as Cardassion 
term by the authors) may arise as a consequence of embedding our observable universe as a 
(3+l)-dimensional brane in extra dimensions, though an elegant and unique 5-dimensional 
energy-momentum tensor T^j^, that gives rise to equation (1) has not yet been found. If 
n < 1, the new term dominates at late times which implies a modification of gravity at very 
low energy scales. Particularly, if n < 2/3, it gives rise to a positive acceleration. Note 
that, in this scenario, although the universe is flat and accelerating, it contains matter (and 
radiation) only without any dark energy contribution. 

The main goal of this paper is to set observational constraints on this Cardassian expan- 
sion model and check whether it is consistent with current cosmological data. We perform 
a combined analysis of data including the dimensionless coordinate distances to type la su- 
pernovae (SNela) and Fanaroff-Riley type lib (FRIIb) radio galaxies compiled by Daly and 
Djorgovski (2003) and the X-ray gas mass fraction in clusters of galaxies published by Allen 
et al. (2002,2003). These results are timely and complementary to the previous constraints 
from the angular size of high- 2; compact radio sources (Zhu and Fujimoto 2002), CMBR (Sen 
and Sen 2003a,b), the SNela database (Zhu and Fujimoto 2003a; Wang et al. 2003; Cao 
2003; Szydlowski and Czaja 2003; Godlowski and Szydlowski 2003), large scale structures 
(Multamaki et al. 2003) and from optical gravitational lensing surveys (Dev, Alcaniz and 
Jain 2003). 

The plan of the paper is as follows. In the next section, we provide a brief summary 
of the Cardassian expansion scenario relevant to our work. Constraints from dimensionless 
coordinate distance data of SNela and FRIIb radio galaxies are discussed in section 3. In 
section 4 we discuss the bounds imposed by X-ray gas mass fraction in galaxy clusters. 
Finally, we present a combined analysis, our concluding remarks and discussion in section 5. 



In the modified Friedmann equation Eq.(l), there are two model parameters, the power 
index n and the coefficient C of the Cardassian term. Instead of C, it is convenient to use 
the redshift Zcax, at which the two terms inside the bracket of Eq.(l) are equal, as the second 
parameter of the model, i.e, C = [p(-Zcar)]^~"- If we further ignore the radiation components 
in the universe that are not important for the cosmological tests considered in this work, we 
have p = po(l + z)^, p(-Zcar) = Po(l + -^car)^ and 



2. The Ceirdassian model: basic equations 
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where po is the present matter density of the universe. Hence Evaluating the Hubble pa- 
rameter today gives = ^y^pofl + (1 + -^car)^''^""'^]- Conventionally, the critical density of 
the universe is pc = SHq/ShG = 1.8791 x 10^^^/i^g cm~^, where h is the present day Hubble 
constant in units of 100 km s^^Mpc^^, and the present matter density of the universe is 
written in terms of pc as po = ^mPc where Qm is the standard matter density parameter. 
Now in Cardassian model, matter alone makes the universe flat, which means that po = Pc,car 
and Pc,car is the critical density of the universe in Cardassian expansion model (Preese and 
Lewis 2002) 

^ ^'^''^^^ ^ 87rG[l + (1 + ^car)3(^-")] ■ 

and the standard matter density parameter of the universe is — [l+(l+-2^car)^^^~"'^]~^ (note 
that ilm,car = Po/Pc,car = 1, as cxpcctcd for a flat universe). As it shown in equation (3), for 
some combinations of the parameters n and Zcm-, the critical density of the Cardassion model 
can be much lower than the one of the standard Fricdmann model. In other words, in the 
context of Cardassian model, it is possible to make the dynamical estimates of the quantity 
of matter that consistently point to po — (0.2 — 0.4)pc compatible with the observational 
evidence for a flat universe from CMB observations and the flatness prediction made by 
inflationary cosmology without any dark energy component (see Preese and Lewis 2002 for 
more details). 

Now we evaluate the dimensionless coordinate distance, y{z), the angular diameter 
distance, D^{z), and the luminosity distance, D^{z), as a function of redshift z as well as 
the parameters of the Cardassian model. The three distances are simply related to each 
other by = {l + zfD^ = {c/ Ho){l + z)y{z). Following the notation of Peebles (1993), we 
define the redshift dependence of the Hubble paramter H as H{z) — HqE{z). Parametrizing 
the model as {^m-i n), we get E function as 

E\z; Q^, n) = Q^il + zf + (1 - + zf^ (4) 

We note that the quintessence models with a constant equation of state (p = cup) for the 
dark energy component give rise to the same H{z) presented here. One can make the 
following identification: uj = n ~ 1, such that n = 1 corrcponds to a ACDM model. As far 
as any observation that involves only H{z), the two models predict the same effects on the 
observation (Preese 2003). It is straightforward to show that the distances are given by 

c . r dz' 



D''{z;Ho,nm,n) = {l+zf-D\z;Ho,nm,n) = -^{l+z)-y{z-^l^,n) = ^{l+z)- f 

-no -no Jo 



E(z'\Q.jn.n) 
(5) ■ 
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3. Constraints from the dimensionless coordinate distance data 

Recently, Daly and Djorgovski (2003) compiled a large database of the dimensionless 
coordinate distance measurements estimated from the observations of SNela and FR lib 
radio galaxies. The author used the database to derive the expansion rate of the universe as 
a function of redshift, E{z), and the acceleration rate of the universe as a function of redshift, 
q{z) (Daly and Djorgovski 2003). We use this sample to give an observational constraint on 
the model parameters, n and ^rn- 

The SNela measurements include the 54 supernovae in the "primary fit C" used by 
Perlmutter et al. (1999), the 37 supernovae published by Riess et al. (1998), and the 
so far highest redshift supernova 1997ff presented by Reiss et al. (2001). The standard 
procedures of Perlmutter et al. (1999) and Riess et al. (1998) were used to determine the 

dimensionless coordinate distances to the supernovae. The apparent bolometric magnitude 
■m{z) of a standard candle with absolute bolometric magnitude M is related to the luminosity 
distance by m = M + 5 logD^ + 25. Then using the relation of equation (5) the B-band 
magnitude-redshift relation can be written as 

TUB ^ Mb + 5\og[c(l + z) ■ y(z)] (6) 

where Ais = Mb — 5 log Hq + 25 is the "Hubble-constant-free" 5-band absolute magnitude 
at maximum of a type la supernova (SNIa). Then the above relation is used to determine 
y{z) for each SNIa. There are 14 SNela that are present in both the Perlmutter et al. (1999) 
and Riess et al. (1998) samples, for which we will use the average values of y with appropriate 
error bars (see Table 4 of Daly and Djorgovski 2003). Therefore we totally have 78 SNIa 
data points that are shown as sohd circles in Figures 1 and 2. 

The use of FRIIb radio galaxies to determine the angular size distance or dimensionless 
coordinate distance of them at different redshifts was first proposed by Daly (1994) (see 
also Guerra, Daly, and Wan 2000; Daly and Guerra 2002; Podariu et al. 2003; Daly and 
Djorgovski 2003). This method consists in a comparison of two independent measurements 
of the average size of the lobe-lobe separation of FR lib sources, namely, the mean size 
< D > of the full population of radio galaxies at similar redshift and the source average 
(over its entire life) size D^,, that is determined via a physical model describing the evolution 
of the sources. The basic idea is that < D > must track the value of -D*, such that the ratio 
i?* =< D > /D^ is assumed to be a constant: R^,{P,y{z)) = k. ., where P is one parameter 
entering into the ratio i?*. y{z) can be determined using an iterative technique, as described 
in detail by Guerra, Daly and Wan (2000) and Daly and Guerra (2002). We use the values 
of y{z) for 20 FRIIb radio galaxies obtained using the best fit to both the radio galaxy and 
supernova data (see Table 1 of Daly and Djorgovski 2003), that are shown as empty squares 
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in Figure 1 and 2. The best fit values of k and (3 and their error bars are included in Table 
2 of Daly and Djorgovski (2003), i.e., k = 8.81 ± 0.05 and /3 = 1.75 ± 0.04. In determining 
the error bar on y{z), the uncertainties of k, and f3 have been included (Daly and Djorgovski 
2003). It is very important to consider whether significant covariance exists between the 
different parameters determined by the fit (Daly and Guerra 2002). In the framework of 
quintessence model, Daly and Guerra (2002) have estimated the likelihood contours in the 
(5 — uj plane and the (5 — Vtm plane respectively (see Figure 3 and 4 of their paper). As 
was pointed out in section 2, the Cardassion scenario is equivalent to quintessence model by 
identifying uj — n — 1. Therefore their results are exactly appropriate to our parameters n 
and Q^: i.e., there is no covariance between (5 and n (^1^). 

We determine the model parameters n and Vtrn through a minimization method. The 
range of n spans the interval [-1,1] in steps of 0.01, while the range of 0,^ spans the interval 
[0, 1] also in steps of 0.01. 



where y{zi]n,Vtrn) refers to the theoretical prediction from equation (5), yoi is the observed 
dimensionless coordinate distances of SNela and FRIIb radio galaxies, and CTj is the uncer- 
tainty (i refers to the ith data point, with totally 98 data). The summation is over all of the 
observational data points. 

The results of our analysis for the Cardassian expansion model are displayed in Figure 
3. We show 68.3% and 95.4% confidence level contours in the (r2m,n) plane using the 
lower shaded and the lower plus darker shaded areas respectively. The best fit happans at 
Q.jn — 0.38 and n — —0.20. It is clear from the figure, that the dimensionless coordinate 
distance test alone weakly constrains the Cardassian expansion model. Only models with 
VLm > 0.60 and n > 0.54 are excluded at the 95.4% confidence level. However, this already 
strongly suggest an accelerating universe (because n < 0.54 < 2/3). Moreover, as we shall 
see in Sec. 5, when we combine this test with the X-ray gas mass fraction test, we could get 
very stringent constraints on the Cardassian scenario. 

Recently, Knop et al. (2003) investigated in detail the effects of various systematic errors 
of SNela on the cosmological measurements (section 5 of their paper). Their main results 
were summarized in Table 9 (Knop et al. 2003), for example, the differences in lightcurve 
fitting methods can change the fiat universe value of flm by ~ 0.03 and o; (= n — 1) by 
0.02. Other systematic errors considered include non-type la SN contamination, Malmquist 
bias, X-corrections, SN colors, dust extinction and gravitational lensing etc. All identified 
systematic errors together give rise to AQ^ = 0.04 and Auj {— An) = 0.09, which are smaller 
than than the current statistical uncertainties of SNela (Knop et al. 2003). 




(7) 
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4. Constraints from the galatxy clusters X-ray data 

Clusters of galaxies are the largest virialized systems in the universe, and their masses 
can be estimated by X-ray and optical observations, as well as gravitational lensing mea- 
surements. A comparison of the gas mass fraction, /gas = M^g^/Mi^t, as inferred from 
X-ray observations of clusters of galaxies to the cosmic baryon fraction can provide a direct 
constraint on the density parameter of the universe Qrn (White et. al. 1993). Moreover, 
assuming the gas mass fraction is constant in cosmic time, Sasaki (1996) show that the /gas 
measurements of clusters of galaxies at different redshifts also provide an efficient way to 
constrain other cosmological parameters decribing the geometry of the universe. This is 
based on the fact that the measured /gas values for each cluster of galaxies depend on the 
assumed angular diameter distances to the sources as /gas oc [D'^]^/^. The ture, underlying 
cosmology should be the one which make these measured /gas values to be invariant with 
redshift (Sasaki 1996; Allen at al. 2003). 

Using the Chandra observational data, Allen et al. (2002; 2003) have got the /gas profiles 
for the 10 relaxed clusters. Except for Abell 963, the /gas profiles of the other 9 clusters appear 
to have converged or be close to converging with a canonical radius r25oo, which is defined 
as the radius within which the mean mass density is 2500 times the critical density of the 
universe at the redshift of the cluster (Allen et al. 2002, 2003). The gas mass fraction values 
of these nine clusters at r25oo (or at the outermost radii studied for PKS0745-191 and Abell 
478) are shown in Figure 4. We will use this database to constrain the Cardassian expansion 
models. Following Allen et al. (2002), we have the model function as 



/gas i^i'i ^1 ^m) 



60h 



h 



1 3/2 



0.5D4(z,;n,l]^), 



(8) 



where the bias factor b 0.93 (Bialek et al. 2001; Allen et al. 2003) is a parameter motivated 
by gas dynamical simulations, which suggest that the baryon fraction in clusters is slightly 
depressed with respect to the Universe as a whole (Cen and Ostriker 1994; Eke, Navarro 
and Frenk 1998; Frenk et al. 1999; Bialek et al. 2001). The term {h/O.bf/^ represents the 
change in the Hubble parameter from the defaut value of Hq — 50kms~^Mpc~^ and the 
ratio D^Qj^j^{zi)/D^g^^{zi;n,^m) accounts for the deviations of the Cardassian model from 
the default standard cold dark matter (SCDM) cosmology. 

Again, we determine the Cardassian model parameters n and through a min- 
imization method. We constrain n^h^ = 0.0205 ± 0.0018 (O'Meara et al. 2001) and 
h = 0.72 ± 0.08, the final result from the Hubble Key Project by Freedman et al. (2001). 
The range of n spans the interval [-1,1] in steps of 0.01, while the range of fi^ spans the 
interval [0, 1] also in steps of 0.01. The difference between the model function and SCDM 
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data is then (Allen et al. 2003) 



i=l 




+ 



h - 0.72 
0.08 



(9) 



where f^°'^{zi;n,flm) refers to equation (8), /gas.oi is the measured /gas with the defaut 
cosmology (SCDM and Hq = 50km Mpc~^) and 0"/^^^ - is the symmetric root-mean-square 
errors {i refers to the ith data point, with totally 9 data). The summation is over all of the 
observational data points. 

The results of our analysis for the Cardassian expansion model arc displayed in Figure 
5. We show 68.3% and 95.4% confidence level contours in the {Qrn,n) plane using the 
lower shaded and the lower phis darker shaded areas respectively. The best fit happans at 
Qm = 0.30 and n = 0.14. It is clear from the figure, that although the X-ray gas mass 
fraction test alone constrains the density parameter very stringently, it weakly limits the 
Cardassian power index n. However, when comparing Figure 4 with Figure 3, we can expect 
the X-ray gas mass fraction test to break the degeneracy presented in the dimensionless 
coordinate distance test of last section. 

As Figure 5 shown, measurements of /gas of galaxy clusters provide an efficient way to 
determine Q^. However the uncertainty of the bias factor b will lead to a systematic error in 
this kind of analysis (Allen et al. 2003). Because it hnearly scales the X-ray mass fraction 
/gas in equation (8), lowering (uppering) it by ~ 10% would cause the best fitting value of 
flm to reduce (increase) by a similar amount. Another systematic uncertainty comes from 
the /gas profiles of galaxy clusters: any rise in the /gas values beyond the measuremet radii 
would cause a corresponding reduction in flm. 



Now we present our combined analysis of the constraints discussed in the previous sec- 
tions and summarize our results. In Figure 6, we display the 68.3%, 95.4% and 99.7% 
confidence level contours in the {Qrn,n) plane from a combination of databases of dimen- 
sionless coordinate distances to SNcla and FRllb radio galaxies and the X-ray gas mass 
fraction in clusters of galaxies. The best fit happans at Qm = 0.30 and n = 0.06. As it 
shown, although the two Cardassion parameters are not very sensitive to the dimensionless 
coordinate distance data of SNela and FRIIb radio galaxies and the Cardassion power index 
is not sensitive to the X-ray gas mass fraction data of clusters, a combination of the two 
data sets gives at the 95.4% confidence level that flrn — 0.301q;q2 and n — O.Oeloiig; ^ very 
stringent constraint on the Cardassion expansion scenario. 



5. Combined analysis, concluding remarks and discussion 
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Given the two model parameters, n, and Q^, Zhu and Fujimoto (2003b) derived the 
redshift Zcax and Zg=o satisfying the relation as follows 



(1 + z\=o = (2 - 3n)^{l + z,^) 



(2-3n)(-^-l) 

u uiyy-i 



3(l-7l) 

(10) 



where z^ar is the redshift at which the two terms inside the bracket of Eq.(l) arc equal, while 
Zq=Q is the redshift at which the universe switches from deceleration to acceleration, or in 
other words the redshift at which the deceleration parameter vanishes. It was shown (Zhu 
and Fujimoto 2003b) that, for every value of VLm, there exists a value for the power index 
of the Cardassian term, n^eski^m), satisfying (2 - 3npeak)~^ exp [3(1 - n^esk)/ (2 - 3npeak)] = 
— 1 which makes the turnaround redshift Zq=Q reach the maximum value, [2:5=0] max = 
exp [1/(2 — 3npeak)] — 1- Lower is, higher [2;g=o]max will be. For the lower bound obtained 
here, Jl^ = 0.28, we have ripeak = 0.0576 and [-2y=o]max = 0.73. In conclusion, our combined 
analysis result, Vtm = 0.30lgo2 ^"^^ ^ — ^-^^^oil 95.4% confidence level, suggest 

a Cardassion expansion universe which swithes from deceleration to acceleration around 
Zq=o e (0.52, 0.73). However, the modified term of the Friedmann equation would dominate 
at a redshift around Zca,v ^ (0.25,0.55), a little bit later than the expansion turnaround 
happans (note that, ^car is generally not equal to Zg=o), which is simply due to the resulting 
power index n is well below 1/3 (Zhu and Fujimoto 2003b). 

Alternative cosmologies from a modified Friedmann equation may provide a possible 
mechanism for the present acceleration of the universe congruously suggested by various 
cosmological observations. In this paper we have focused our attention on one of these 
scenarios, the so-called Cardassion expansion in which the universe is flat, matter (and 
radiation) dominated and accelerating but without any dark energy component. We have 
shown that stringcnct constraints on the parameters n and flm, that completely characterize 
the scenario, can be obtained from the combination analysis of the dimensionless coordinate 
distance data of SNela and FRIIb radio galaxies and the X-ray mass fraction data of clusters. 
It is natually hopeful that, with a more general analysis such as a joint investigation on 
various cosmological observations, one could show clearly if this scenario constitutes a feasible 
alternative to other acceleration mechanisms. 
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Fig. 1. — Dimensionlcss coordinate distances y{z) as a function of log 2; for 78 SNela and 20 
FRIIb radio galaxies. The solid circles mark the SNela, while the empty squares mark the 
FRIIb radio galaxies. The solid curve corresponds to our best fit to the total 98 data points 
with — 0.38, n = —0.20. The database are taken from Daly and Djorgovski (2003). 
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Fig. 2. — The residuals of y{z) relative to our best fit model with = 0.38, and n — —0.20 
as a function of log z for 78 SNela (the sohd circles) and 20 FRIIb radio galaxies (the empty 
squares). The data are taken from Daly and Djorgovski (2003). 
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Fig. 3. — Confidence region plot of the best fit to the database of the dimensionless coordinate 
distances to 78 SNela and 20 FRIIb radio galaxies compiled by Daly and Djorgovski (2003) - 
see the text for a detailed description of the method. The 68.3% and 95.4% confidence levels 
in the n-Q^ plane are shown in lower shaded and lower + darker shaded areas respectively. 
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Fig. 4. — The apparent redshift dependence of the /gas measured at r25oo for 9 clusters of 
galaxies with convergent /gas profiles. The error bars are the symmetric root-mean-square 
1(7 errors. The solid circles mark the six clusters studied by Allen et al.(2002), while the 
empty squares mark the other three clusters published by Allen et al. (2003). The solid 
curve corresponds our best fit to the Cardassian model with flrn = 0.30, and n = 0.14. 
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Fig. 5. — Confidence region plot of the best fit to the /gas of 9 clusters published by Allen 

et al. (2002,2003) - see the text for a detailed description of the method. The 68.3% and 
95.4% confidence levels in the n-flm plane are shown in lower shaded and lower + darker 
shaded areas respectively. 
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Fig. 6. — Confidence region plot of tfie best fit from a combined analysis for the dimension- 
less coordinate distances to 78 SNela and 20 FRIIb radio galaxies compiled by Daly and 
Djorgovski (2003) and the X-ray gas mass fractions of 9 clusters published by Allen et al. 
(2002, 2003). The 68.3%, 95.4% and 99.7% confidence levels in the n-Q^ plane are shown 
in white, white + lower shaded and white + lower and darker shaded areas respectively. 



